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Abstract. We present preliminary results from a series of high-resolution 
A^-body simulations that focus on 8 dark matter halos each of order a 
million particles within the virial radius. We follow the time evolution of 
hundreds of individually tracked satellite galaxies and relate their phys- 
ical properties to the differing halo environmental conditions. Our main 
science driver is to understand how satellite galaxies lose their mass and 
react to tidal stripping. Unlike previous work our results are performed 
in a fully self-consistent cosmological context. The preliminary results 
demonstrate that while environment may vary considerably with respects 
to formation time and richness of substructure, the satellites evolve sim- 
ilarly. 



1. The Simulations 

Four standard ACDM cosmological simulations with side length 64/i~^ Mpc were 
run at low mass resolution (128'^ particles), with 8 halos selected sampling a 
range of environments. In this contribution we present results for our oldest 
(8.3 Gyrs) and youngest (3.4 Gyrs) dark matter halos, halo # 1 and # 8 respec- 
tively. These halos were re-simulated at 64 times higher mass resolution giving 
an effective mass per particle of 10^ Mq, and a force resolution of approxi- 
mately 1 kpc. The A^-body simulations were carried out using the publicly 
available adaptive mesh refinement code MLAPM (Knebe, Green & Binney 2001). 
The code solves Poisson's equation on a hierarchy of nested grids: the whole 
computational volume is covered by one cubic domain grid whereas refined re- 
gions are of arbitrary shape and adjusted to the actual density field at each 
major step in order to follow the real distribution of particles at all times. An 
example of MLAPM in action is shown in the left panel of Figure 1. These two im- 
ages show the density field of our youngest halo superimposed by the refinement 
grids. The top image shows the 5th refinement level with grid resolution 4096^, 
the bottom the 6th level of refinement. As the grids follow the over-densities in 
the simulations they by definition encompass the satellites. Hence, we can build 
a hierarchy of nested isolated MLAPM grids into a " grid tree" , and then generate 
a list of potential satellite centers by treating the densest grid at the end of each 
branch as a "potential" halo center. We further check the validity of the center 
by assigning physical properties to the adjacent radial density distribution. The 
details can be found elsewhere (Gill, Knebe & Gibson 2003). 
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Figure 1. (Left Panel) Consecutive refinement levels of MLAPM's grid structure 
for halo # 8 superimposed upon the density projection of the particle distribution. 
(Right Panel) Orbits of all objects in the vicinity of the host halo; orbital lines 
graduate from dark at Zform time to light z = 0. 



2. Dark Matter Halo Environment 

Following Lacey & Cole (1994) the formation time of the halo Zform is defined 
to be the time when the halo contains half of its present day mass. From ^form 
to z=0 we track all satellites in the vicinity of the host halo. The number of 
satellites within the virial radius of our 8 halos (of mass > 2 x 10-'^'^ Mq) 
ranges from 50 — 250. The right panel of Figure 1 shows the orbits of the tracked 
satellites for two sample host halos highlighting the differing halo environments. 
The top image shows a halo developing in a rich region fed by three filaments, 
where as the figure directly below shows a host halo in a relatively isolated region 
that saw early rapid collapse. 



3. Satellite Disruption 

Figure 2 shows the number of satellites within a 2 Mpc sphere (physical 
coordinates), as a function of redshift normalized by the number of satellites 
within the sphere at redshift ;Zform- The thin line (always on top) shows the 
total number of satellites within that sphere, while the bottom (thick) line shows 
the number of "disrupted" satellites subtracted from the number of satellites 
at Zform- The rise in the thin lines from Zformto z=0 reflects the richness of 
the immediate environment around the halo. Halo #1 is in a particularly rich 
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Figure 2. (Left Panel) Number of satellites within a 2h~^ Mpc real co-ordinate 
sphere, as a ftinction of redshift. Please note that the a;-axis' are not the same range 
but rather covering the time from redshift ^; = to 2:forin for each individual halo. 
The thin line is the total number of satellites within that sphere whereas the thick 
line describes the disruption rate of the satellites within the sphere. (Right Panel) 
Mass loss per orbit for satellites that survived until redshift 2: = as a function of 
their orbital eccentricity e. 



"filament fed region". These filaments are a source of significant satellite infall 
(recall the right panel of Figure 1). 

As the satellites orbit within the halo they experience dynamical friction 
causing them to sink towards the centre of the halo's potential, undergo tidal 
stripping, and gradually lose mass. For each satellite, its particle distribution 
becomes increasingly diffuse. Our criterion for calling a satellite "disrupted" is 
linked to its tidal radius; assuming that the average density of the satellite has 
to be at least a factor of three above the density of the host halo at the distance 
of the satellite we iteratively calculate the tidal radius and the mass m(< rtidai) 
enclosed by solving: 

rt.dai - ( 3^(^ ^) ) D = 0, (1) 

The distance D of the satellite to the host and M(< D), the mass of the host 
halo interior to D, arc fixed. As soon as there are fewer than 15 particles within 
the tidal radius we classify the satellite as being disrupted. Note that we are 
unable to separate numerical resolution disruption and real physical disruption 
of a satellite. It is not clear if we had infinite mass resolution that the satellite 
would still actually survive. 

In Figure 2, the right panel shows the difference of the satellite's mass at 
Zformand z=0 divided by its mass at Zformand the number of orbits as a function 
of eccentricity of the satellite orbit. Thus as the satellite does get disrupted in 
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general it losses at least 30 percent of its mass each orbit. Further, for seven of 
the eight halos, satellites on high eccentricity orbits are losing more mass (per 
orbit) than spherical ones, (such as halo # 8). However, this is not strictly true 
for all environments; e.g. in halo # 1, satellites on spherical orbits lose more 
mass - in this case the satellite orbits are on average closer to the halo core. 



4. Properties of the Satellites in the Halos 

Because of the high temporal resolution of the outputs of our simulations (At ~ 
0.2 Gyr) we are able to track accurately the orbital and other characteristics of 
the satellites. We present here a few preliminary results: 

• For each of the 8 halos the mean eccentricity was consistently around 0.6, 
where eccentricity was measured as e = 1 — ap"cc°tre • 

• For each of the 8 halos the mean pericentre distance was approximately 
30 percent of the host virial radius. 

• The apocentres of the satellite orbits are aligned with the principle axis of 
the host galaxy cluster. 

• In general a satellite loses at least 30 percent of its mass each orbit. Fur- 
ther, for the majority of the eight halos, satellites on high eccentricity 
orbits are losing more mass (per orbit) because they penetrate deeper into 
the host's potential. 

• Satellite-Satellite interactions are not negligible. 

A fuller and quantitative analysis is being prepared (Gill, Knebe & Gibson 
2003). There, we address additional characteristics concerning satellite evolu- 
tion: e.g. how does a satellite actually react to its tidal disruption? 
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